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Abstract 

Unitary evolution from pure initial states to pure final states in 7r~p — » 7r~7r"'"n imposes constraints on 
pion production amplitudes that are violated by the CERN data on polarized target at 17.2 GeV/c (isj . 
The pion creation process is a non- unitary evolution of initial state pi{n~p) to final state p/(7r^7r+n) arising 
from a unitary co-evolution of the pion creation process with a quantum environment described by Kraus 
representation. The purpose of this work is to identify the interacting degrees of freedom of the environment 
in a high resolution amplitude analysis of CERN data on n^p Tr^vr+n on polarized target for dipion masses 
580-1080 MeV where S- and P-wave transversity amplitudes dominate. The 5'-wave spectra show presence of 
p'~'{770) arising from the presence of p°(770) peak in data component ai + a2 in relation |5p = oi -|-a2 — 3|Lp 
while the P-wave spectra |Lp show a dip at /o(980) mass arising from /o(980) structures at 980 MeV in 
all data. The observed p"(770) — /o(980) mixing is encoded in all measured density matrix elements which 
also encode a level splitting of the spectra arising from the interaction of the pion creation process with the 
environment. The analytical form of the level splitting reveals the existence of a new quantum number g 
characterizing the environment and allows to identify the four sets of solutions for the transversity amplitudes 
with the four co-evolution amplitudes required by the Kraus representation. We propose a model for the 
CPT violating and non-dissipative interaction of the pion creation process with the environment in which 
non-diagonal transitions between resonant qq modes and 7r"7r+ states lead to vector-scalar mixing. The 
model predicts quantum entanglement of 7r~7r"'" isospin states. The final states in 7r~j5 — > 7r~7r"'"n do not 
posses prepareable CPT conjugate states due to entanglement of tt^tt^ pairs and because the environment 
states do not have well defined charge conjugate states. As a result the concept of CPT symmetry looses 
its meaning in the unitary co-evolution of the pion creation process with quantum environment. 
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I. INTRODUCTION - THE PUZZLE OF cr(750) RESONANCE. 



Following the discovery in 1961 of p meson in vrA^ — > ttt^N reactions (Ij, the measurements of 
forward-backward asymmetry in tt~p — > Tr^vr+n suggested the existence of a rho-like resonance in 
the S'-wave, later referred to as (t(750) scalar meson [2, H, 13) @] • It was expected that the (t(750) 
resonance would show up prominently in 7r~p — > Tr^n^n production where the ^^^(770) does not 
contribute. However the measurements of this reaction at CERN in 1972 found no evidence for 
a rho-like (t(750) [7]. Furthermore, in 1973, Pennington and Protopopescu used analyticity and 
unitarity constraints on partial wave amplitudes in tttt vrvr scattering (Roy equations) to show 
that a narrow (t(750) resonance cannot contribute to vrvr scattering [8]. From these facts it was 
concluded that (t(750) does not exist and in 1974 Particle Data Group dropped this state from its 
listings. 

In 1972, van Rossum and his spin physics group at Saclay reported the first measurements of 
recoil nucleon polarization in vrA^ — > vrA^ elastic scattering made at CERN at 6 and 16 GeV/c 0]. 
The resulting complete set of observables enabled the first model independent amplitude analysis 
of a hadronic reaction [l^] . The results invalidated all Regge models (but not the concept of Regge 
poles). These findings established that the experimental determination of amplitudes in hadron 
scattering using measurements with spin is crucial for our understanding of hadron dynamics. 



In 1978, Lutz and Rybicki extended the concept of amplitude analysis to pion production 
processes. They showed ll| that almost complete amplitude analysis of reactions vrA^ tt'^tt~N 
and KN K^tt~N is possible from measurements in a single experiment on a transversely 
polarized target. The work of Lutz and Rybicki opened a new approach to hadron spectroscopy 
and hadron dynamics by enabling us to study the production of resonances on the level of spin 
amplitudes rather than spin- averaged cross-sections. 



The pion production on polarized targets was measured at CERN in several processes. 
CERN-Munich-Cracow group measured 7r~p — > TT~TT'^n on transversely polarized proton target 
in a high statistics experiment at 17.2 GeV/c at low momentum transfers —t = 0.005 — 0.20 
{GeV/cf [l3, [H, Ei] and at high momentum transfers -t = 0.20 - 1.00 (GeV/c)^ Saclay 
group measured n^n vr+vr^p 0, [13] and K+n K+ir-p 0, [ij on transversely polarized 
deuteron target at 5.98 and 11.85 GeV/c at larger momentum transfers —t = 0.20 — 0.40 (GeV/c)^. 
Recently measurements of vr^p 7r~7r^n on transversely polarized proton target at 1.78 GeV/c 
at low momentum transfers —t = 0.005 — 0.20 (GeV/c)^ were made at ITEP [201] . 



The CERN measurements of vr p ^ it vr+n and vr+n tt^tt p on polarized targets reopened 
the question of the existence of (t(750) scalar meson. Evidence for a narrow (7(750) was found 



in amplitude analyses of n p 



m, Ha, 12J, \M and in vr 



n 



vr+TT p at 



vr-vr+n at 17.2 GeV/c 

5.98 and 11.85 [25|,l2d]. Clear evidence for o-(750) emerged from later and more precise amplitude 
analyses of both reactions 27, 2^, 2^. New evidence for (t(750) comes from the amplitude analysis 
of the ITEP data at 1.78 GeV/c [3. The best latest fit in Ref. [29( gives = 778 ± 16 MeV 
and Tu = 142 it 33 MeV. These values are very close to resonance parameters of p(770) with 



nir. 



775.8 ± 0.5 MeV and = 150.3 ± 1.6 MeV. 



In 2001, E852 Collaboration at BNL reported high statistics measurements vr p — > vr'^vr^n at 
18.3 GeV/c [30]. The data revealed large differences in the S'-wave intensities in ■k~-k~^ and 
production. There was no evidence for a rho-like a{750) resonance in the vr'^Tr'' S'-wave. The BNL 
data presented anew the puzzle of the (t(750) resonance. The CERN and BNL data are both 
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high quahty data that cannot be used to exclude one another. We must accept that they are 
both correct and that the apparent contradictions between them are tehing us something new and 
important. The puzzle of cr(750) resonance has become a unique opportunity to learn new physics. 



The first hint on the solution of the puzzle emerged from the new amplitude analysis of CERN 
measurements of 'K~p — > 'K~'n~^n in Ref. 29(]. This study extended the dipion mass range of 600 
- 900 MeV of previous analyses 0, ^ to 580 - 1080 MeV to include the /o(980) resonance. 
Both solutions for the P-wave amplitude \L(i\'^ showed an unexpected dip at ~ 980 MeV, the 
mass of /o(980) resonance. In 2001, Gale and collaborators investigated p(770) — ao(980) mixing 
in lepton pair production at RHIC _arising from the violation of Lorentz symmetry by vacuum 

|2 



31C. 



and found large measurable effects 
p{770) — /o(980) mixing in the P-wave. 



Their work suggested that the dip in arises from 

Similar mixing in the S-wave then allows to identify 
cr(750) with p^{770) which immediately explains why no (t(750) was ever found in ir^p 
production. 



The S-wave intensities in vr+n — > vr+vr^n at 11.85 GeV/c and iT~p tt^tt^ at 18.3 GeV/c show 
particularly large differences in their structure for masses below ~ 600 MeV. The expectation 
that these intensities should be similar in this mass range can be traced to the isospin structure 
of the amplitudes which originates in the assumption of Generalized Bose-Einstein symmetry 
for two-pion states. The assumption of Generalized Bose-Einstein symmetry leads to three 
relations that must be satisfied by combinations of partial wave intensities of 7r~7r+, vr'^vr'^ and 
ir~^7r~^ production for even dipion spins. In a sequel paper we will show that the relations are all 
violated by the available data indicating a violation of Generalized Bose-Einstein symmetry in 

7r~p TT~TT~^n. 

There are four sets of physical solutions for the S- and P- wave transversity amplitudes Au{i) and 
^d(j)) ^) J = 1) 2 with target nucleon transversity "up" (u) and "down" (d) relative to the scattering 
plane. All these solutions involve p''(770) — /o(980) mixing and encode relevant information about 
the pion creation dynamics. The multitude of the physical solutions may reflect important new 
aspects of reality. To select one of these solutions as the only valid physical solution amounts to a 
loss of relevant information. We can retain the full information about the dynamics encoded in the 
amplitudes from all solutions if we accept that all these solutions are valid physical solutions. But 
then the measured final state density matrix pj must be a mixed state formed by density matrices 
Pf{ij) with probabilities pij corresponding to all four solutions i,j = 1,2 

Pf{Tr~Tr+n) = pnpfill) + Pi2pf{l2) + p2ipfi2l) +p22Pf{22) (1.1) 

where ^Pij = 1- The final state is a mixed state even when all states Pf{ij) are pure states 
arising from pure initial states. If the hypothesis (1.1) is true then pure initial states would evolve 
into mixed final states in TT~p iT~Ti~^n and the dynamics of the pion creation process would be 
a non-unitary evolution. The non-unitary dynamics could be responsible for the p^{77Q) — /o(980) 
mixing and the violations of Generalized Bose-Einstein symmetry. 



In 1980 Wald showed that any scattering process of particles that evolves pure initial state into 
a mixed final state violates CPT symmetry and is time irreversible [s^]- In 1982 Hawking showed 
that in the presence of a black hole (macroscopic or microscopic) a pure initial state of interacting 
particles will evolve into a mixed final state as some of the quantum states produced in the 
particle interaction will fall behind the horizon and become inaccessible to measurements by the 
observer [s^] • He also suggested that quantum fluctuations of the space-time metric will have the 



same eff'ect on interacting particles and induce their non- unitary evolution - at any energy 33|, |3j] . 
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Hawking questioned the universal validity of the unitary time evolution in Quantum Field Theory 
in the presence of metric fluctuations and suggested that initial and final state density matrices 
Pin and Pout are connected by a linear but non-unitary evolution superoperator. 

Hawking's ideas inspired suggestions to test them experimentally. In 1984, Ellis, Hagelin, 

Nanopoulos and Srednicki proposed that quantum fluctuations of the metric form an environment 

with which interacting as well as free particles interact as open quantum systems [35]. Such interac- 

n — 

tions would lead to an observable change of coherence and CPT violations in K K systems. Over 
the years other suggestions to test Hawking's ideas have been put forward 36l.l37l.l3 



Recent experiments with free neutral kaons have yielded remarkably sensitive results on violations 



of CPT symmetry, time reversal invariance and entanglement of kaon pairs [43|, |4J]. So far 
these experiments did not provide a conclusive confirmation of a non-unitary evolution, and thus 
possible evidence for quantum gravity effects. 



In our previous work [45[ we returned to the original Hawking's idea and searched for evidence 
of a non-unitary evolution of pure initial states of interacting particles into mixed final states in 
the existing CERN data on vr^p vr^vr^n on polarized target at 17.2 GeV/c. The purity of the 
final state pf is controlled by the purity of the recoil nucleon polarization. We developed a spin 
formalism to calculate the expressions for recoil nucleon polarization for two specific measured 
initial pure states. Imposing the condition of purity on the recoil nucleon polarization we obtained 
conditions on the amplitudes which are violated by the model independent amplitude analyses of 
the CERN data on polarized target at large momentum transfers. We conclude that pure states 
can evolve into mixed states in vr^p — > 7r~7r"'"n. 

In quantum theory such non-unitary evolution occurs in open quantum systems or processes S 



interacting with a quantum environment E [47|, |48|, |49| . The co-evolution of the system S and the 
environment E is unitary 

Pf{S, E) = Upi{S) ® pi{E)U+ (1.2) 

The reduced density matrix 

Pf{S)=TrE{pf{S,E))=£{p,{S)) (1.3) 

is a mixed state and £ describes the non-unitary evolution from pi{S) to Pf{S). It is given by 
Kraus representation 47, 4^, 4^ 



Pf{S) =£{pi{S)) = J]]^Pmn5'£mPi(5')S'+^ (1.4) 



i m,n 



where Sim =< G£\U\em > and |e^ > are interacting degrees of freedom of the environment. The 
initial state of the environment has a general form 

Pi{E) = '^Pmn\em >< e„| (1.5) 

m,n 

We showed in [is^ that Kraus representation leaves invariant the formalism used in data analyses 
provided the co-evolution with the environment conserves P-parity and quantum numbers of the 



environment. On general grounds [45|, |43| there are four interacting degrees of freedom of the 



environment in vr p ^ tt ir^n. The Kraus representation then has a diagonal form 

4 4 

Pf(^) = ^PuSupi{S)Sj^ = ^pepfii) (1.6) 
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The measured density matrix elements are redefined to be environment-averaged density matrix 
elements. The elements that depend explicitely on the environment are predicted to violate certain 
phase relations. The prediction is in excellent agreement with the CERN data 45|] and validates 
the view of pion creation processes as an open quantum systems interacting with a quantum 
environment. 



The question now arises what are the quantum states \e£ > of the environment. Experimentally, 
the final state in tt~p vr^vr+n must be the mixed state Pf{S) given by (1.6) which has a form 
identical to the mixed state pf{TT~Tr^n) given by (1.1). The Central Hypothesis put forward in 
this work is that the two states are the same states. The hypothesis then alows us to relate the 
solutions for the amplitudes with quantum states of the environment. 

We associate with the two solutions for amplitudes Au{i) and Afi{j) two qubit states \i > and 
|j > where i,j = 1,2 and identify the quantum states \e£ > with two-qubit states \i > \j >. 
Experimentally, the two solutions for S- and P-wave amplitudes Au{i) and Ad{j) originate in 
two solutions for mass spectra |^u(^)P and |^rf(j)p which can be thought of as a form of level 
splitting arising from an interaction with the environment. We show in this work that the level 
splitting is characterized by a new quantum number g = ±1 associated with the environment. 
The two solutions for the amplitudes then reflect their dependence on this quantum number with 
Au{i) = Au{gu) and Ad{j) = Au{gd). The quantum states je^ > are then identified with two-qubit 
states \gugd >■ 

The new quantum number g of the environment gives rise to a new CPT violating interaction 
resposible for the level splitting and the non- unitary evolution in 7r~p TT~TT'^n. We formulate a 
model of such interaction in which the environment states {gudd > interact with a coherent state 
of qq resonant modes produced in an intermediate stage 7r~p qqn . The model explains the 
observed p^{770) — /o(980) mixing in terms of environment induced transitions of qq resonant 
modes with a definite spin K to a superposition of two-pion states with different spins J. The 
model leads naturally to a violation of Generalized Bose-Einstein symmetry and to dynamic 
entanglement of 7r~ it~^ pairs with definite spin J. 

The interaction of pion creation process with the environment can be thought of as a scattering 
of an intermediate coherent state of qq resonant modes carrying energy-momentum and spins 
with particles of the environment carrying quantum entanglement (1.5). There is no exchange of 
energy-momentun between the hadrons and the environment in this non-dissipative, entanglement 
changing interaction. The observed process is time-irreversible and violates CPT symmetry. 

In this work we focus on the S- and P-wave decoherence free subspace of the reduced density 
matrix which is the trace over the unmeasured recoil nucleon spin states. The density matrix 
elements of a decoherence free subspace do not depend explicitely on the interaction with the 
environment. In Section II. we show that all S- and P-wave elements of p^ are expressed in 
terms of reduced transversity amplitudes which do not depend on the relative phase u; between 
the 5- wave transversity amplitudes Su and S^- Phase relations arising from the assumption of 
decoherence free subspace result in cubic equations for P-wave amplitudes l-L^P and which 
render the system analytically solvable from data on matrix elements Rep^ and Repy measured 
with transversely polarized targets. 

In Section III. we present results of high resolution amplitude analysis of CERN data on 
7r~p 7r~TT~^n at 17.2 GeV/c for dipion masses 580-1080 MeV as one of 8 reactions analyzed. 
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The analysis used 1 and 5 million Monte Carlo samplings of the data error volume in each 
mass bin. The analysis with 5 million samplings gives essentially identical results, indicating the 
resolution of mass spectra j^rP)''" = u,d is stable. The relative phases are determined up to a 
sign ambiguity and a relative phase to between 5- wave amplitudes Sd and Su- In Section IV. we 
relate the positivity of density matrix and decoherence assumption to observation that the relative 
phases between Sr and Lr amplitudes do not change sign which allows the unique asignment of 
signs to the relative phases. In Section V. we show how the measurements of matrix elements 
Imp'^ and /mp° resolves the sign ambiguity of the relative phases. In a sequel paper [i^ we show 
that the phase w is uniquely determined by a conversion of transversity amplitudes into helicity 
amplitudes. As a result there are four sets of amplitudes Au{i) and A^ij) consistent with the 
requirement of Kraus representation (1.6). 

In Section VI. we present evidence for ^^^(770) — /o(980) mixing. Resonace peaks of p^{770) 
and /o(980) are clearly resolved in S'-wave amplitudes jS^p in contrast to our previous low 
resolution analyses using 40 000 samplings [ij, 13] • The 5- wave moduli are given by a 
relation |5p = ai + a2 — 3|Lp where ai + 02 is a data component. The presence of ^^^(770) 
in S'-wave spectra \S\'^ arises from the presence of p^{770) peak in the data component oi + 02 
which survives the subtraction of p^{770) peak in the P-wave amplitude 3|Lp. In Section VII. 
we present a test of rotational and Lorentz symmetry in tt~p — > TT~7r~^n. We show that the 
large differences observed in p'^(770) widths in unnatural and natural exchange amplitudes |C/dp 
and \N(i\'^ arise from interference of P-wave amplitudes with helicities A = +1 and —1 and not 
from a violation of Lorentz symmetry - which if true could explain also the p^ (770) — /o (980) mixing. 

In Section VIII. we show that the data exclude the existence of a single physical solution for 
the moduli. Instead the data require that the moduli \At-\'^ have a specific analytical form of two 
roots of the cubic equation which can be interpreted as a level splitting of the spectra. In Section 
IX. we show that the two solutions for the moduli |^rP can be rewritten as interferences of two 
amplitudes ZA,r{gT) with phases Qt^a^t characterized by a new quantum number Qr = ±1. The 
phase \a,t determines the level splitting. The amplitudes Au{i) and ^d(j), i, j = 1,2 correspond 
to amplitudes Au{gu) and ^^(5^) with gu,gd = il- In Section X. we identify the quantum states 
of the environment with two-qubit states \i > \j >, or alternalively, with \gu > \gd > and define 
co-evolution amplitudes of vr^p n^TT'^n process with the environment. 

In Section XI. we formulate a model of interaction of pion creation process with the environment 
to explain p^{770) — /o(980) mixing. In the first stage of pion creation process a coherent state is 
formed of resonant qq modes with spins K{qq) and isospins Ik- The interaction of the coherent 
state with the environment results in transitions of resonant modes with spin K and isospin Ik to 
two-pion states with spins J(7r~7r+) which gives rise to /9o(770) — /o(980) mixing. In Section XII. we 
introduce Generalized Bose-Einstein symmetry conserving and violating co-evolution amplitudes 
corresponding to maximally entangled ■k~'k'^ states. The total co-evolution amplitude is their 
superposition which leads to dynamic entanglement of the observed 7r~7r+ states. The change in 
the entanglement content of the 7r~7r^ states arises entirely from the p''(770) — /o(980) mixing 
which in turn originates in the CPT violating non-local interaction with the environment. The 
paper closes with a brief summary in Section XIII. 
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II. S- AND P-WAVE DECOHERENCE FREE SUBSPACE OF REDUCED DENSITY 
MATRIX IN TT-p 



i^{e^, P) = Tr^=^,{{pf{e<i^, P)ii,) = il{ect>) + pjl{ect>) + Pyi^ym + Pzilm (2.1) 



In previous paper [43] we have developed a general spin formalism for the final state density 
matrix in 7r~p — > 7r~7r+n and similar pion creation processes, and defined nucleon helicity and 
transversity amplitudes with definite t-channel naturality. In the following we assume a familiarity 
with the Section II. of this paper. 

The most feasible experiments are measurements of tt~p — > 7r~7r^n on unpolarized or polarized 
targets with target polarization P = {Px, Py, Pz)- In modern polarized targets the direction of the 
polarization vector P can be selected at will. The experiments measure the two-pion angular distri- 
bution I^{6(j), P) leaving the recoil nucleon polarization vector Q not observed. Such measurements 
provide information on the reduced final state density matrix I^{9(f),P) given by [45| 

i_ 1 

I 2 2 

'xx' 

where the components of I^{6<p,P) are 

72 J max 

^'^^^'^^ = ^ j;^Oj'a(i?ep°)i3f,i?e(y/(e</))y,^*(00)) (2.2) 

J<J'A>0 A' 

for k=u (unpolarized target) and k=y (transversely polarized taget normal to the scattering plane), 
and 

72 J max 

^'^^"^^ = dt£l^ j;Oj'eA(W,)^3f;im(y/(0</.)y/*(00)) (2.3) 

J<J' \>0 A' 

for k=x (transversely polarized target in the scattering plane) and k=z (longitudinally polarized 
target). In (2.2) and (2.3) = 1 and = 2 for A > and the factor ^jj/ = 1 for J = J' and 
^jj' = 2 for J < J' . Experimentally, in a given region of dimeson mass m and momentum transfer 
t only amplitudes with J < Jmax contribute and all sums in (2.2) and (2.3) are finite. It is the 
intensity I^{0(j),P) which has been measured in CERN measurements of pion creation processes 
on transversely polarized targets and from which S- and P-wave density matrix elements were 
determined in small {m,t) bins using maximum likelihood method l3, 16, H3, 51|. 



General expressions for density matrix elements in terms of amplitudes were tabulated by Lutz 
and Rybicki [ll| and are reproduced in Appendix A of Ref. [is! ]. For S- and P-wave amplitudes we 
shall use a simplified notation At- where A = S, L,U, N are transversity amplitudes with definite 
t-channel naturality 

In (2.4) U\ T- ^^'^ ^Xt unnatural and natural exchange transversity amplitudes introduced 
in Ref. [45( and J and A are dimeson spin and helicity. The nucleon transversity t = u,d for 
target spin "up" or "down" relative to the scattering plane, respectively. In Tr~p — > ■K~iT~^n the 
amplitudes C//^ exchange vr and ai quantum numbers in the t-channel while the amplitudes iV/ ^ 
exchange 02 quantum numbers. 

Amplitude analysis of the S- and P-wave subsystem of the reduced density matrix determines 
reduced transversity amplitudes defined as follows 
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TABLE I: Density matrix elements Rep'^cPa / dtdm and RepyCpa/dtdm in terms of reduced transversity 
amplitudes. 

p'l'l' Rep^d^a/dtdm Rep^d^ a / dtdm 

Pfs \S? + \S? 

pII + 

p\\ i(ic/p + ic7n + i(|ivp + \m - + i(|iv|2 _ |iv|2) 
p\u -im'' + \u?) + hm'' + \m -idc^i^ - ic/n + \m'' - \m 

Repl° Re{LS* + LS*) Re{LS* - LS*) 

V2RepYl Re{US* + US*) Re{US* - US*) 

V2Repl\ Re{LU*+LU*) Re{LU* - LU*) 



TABLE IL Density matrix elements Imp'^d^a / dtdm and Imp^d^a/dtdm in terms of reduced transversity 
amplitudes. 



/9^3[ Impld'^a/dtdm Imp^d'^a/dtdm 

V2Impl\ Re{-SN* + NS*) Im(+SN* - NS*) 
V2Impl\ Re{-LN* + NL*) Im{+LN* - NL*) 
Imp]}^^ Re{+UN* - NU*) Im{-UN* + NS*) 



S=\Su\, S=\Sd\ (2.5) 
L = |L„| expi - ^sj , L = \Ld\ expi - ^sj 
U = \Uu\ expi - . U = \Ud\ expi ($[/^ - ^^J 

N = \N^\ expi ($jv„ - , N = lA^dKexp ($jVd - ^sj 

where is the phase of the amplitude At- The reduced transversity amplitudes are related to 
transversity amplitudes by phase factors 

Au = A exp i^Su > Ad = A exp iu exp i^Su (2-6) 

for unnatural exchange amplitudes A = S, L,U and 

Nu = N exp iio exp i^s„ , Nd = N exp i^s„ (2.7) 

for natural exchange amplitude N. In (2.6) and (2.7) is the arbitrary absolute phase and 
uj = ~ ^Su is the relative phase between S-wave amplitudes of opposite transversity. 

In Tables I. and II. we present S- and P-wave density matrix elements expressed in terms of 
the reduced transversity amplitudes. Because of the angular properties of Y^{6(t)), the elements 
(/3^)qq = (p^)s°, (Pfc)oo {pV^vli^ — y independent but appear as two independent 

combinations in the measured angular distributions (2.2) 

{pDsp ^ {pl)f, + {pDH + 2{pl)]i, {pDpp ^ {pDH - (pO)il (2.8) 
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In terms of reduced transversity amplitudes they read 

iPu)sp-^ = \S\' + \S\' + |L|2 + + |[7|2 + |I7|2 + \N\2 + |]V|2 (2.9) 

^ atam 

iPu)pP^ = \L? + \L? - \{\U? + IFP + |iVp + \m (2.10) 
(P°)PP^ = \L?- \L\' - \i\U? + \N\' - \m 

Note that in (2.3) Im{pl)(( = and /m(yo"'((9(/))y/*((90)) = 0. Also note that elements Im{pl)l% 
Im{pl)\^^, Im{pl)\'^, k = u,y and i?e(p^)°];, Re{pl)l\, Re{pl)]}-^^^, k = x,z are not zero. They are 
not observable as the result of parity conservation but they are calculable from amplitude analysis. 

The observables measured in 7r~p — > 7r~7r^n on transversely polarized target organize them- 
selves into two groups involving amplitudes of opposite transversity. Using the expressions in Table 
I. and S = (Pa/dtdm, the two groups are 



ai 



liip'jsp + (/>°)5P)S = \S\' + |L|2 + \U\' + \N\' (2.11) 
a2 = iip'i)pp + ipl)pp)^ = 2\L\'-\U\^-\N\' 
as = iipy.U + ip'y)\'-i)^ = - \U\' 

"4 = l{{p'X, + {pXs)^ = \L\\S\cosi^Ls) 
1 

a, = ^{{p'Xi + ipl)l\)^ = \L\\U\cos{^Lu) 



as = -^{{pl)\', + {p'X)^ = \U\\S\cos{'^us) 



for reduced transversity amplitudes with transversity t = u, and 



ai = 


\{{pI)sp- 




= \s\' + 


|Lp + + \N 


0.2 = 


{{pDpp - 


(pDpp)^ = 




|I7|2 - |ivp 


as = 


{{pl)f-, - 


{P'y)\'-l)^ - 


= \W- 


\W 


04 = 


2((Pu)os ~ 




\L\\S\cos{^Ls) 


as = 


'^^{iPu)ls 


- (pX)^ 


= \U\\S\ 


cos{^us) 


06 = 


"^((Pu)oi 


- ip'y)l\)^ 


= \L\\U\ 


COs($L(7) 



(2.12) 



for reduced transversity amplitudes with transversity t = d. The relative phases are defined as in 

^1) 



(2.5). For dipion masses where 5- and P-wave dominate, {Pu)sp ^-nd {py)sp are traces 



{pl)sp = Tri{p'i)ii) = l, {pl)sp = Tr{{pl)ii)=T (2.13) 
where T is target spin asymmetry [4^. 
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The 6 equations (2.11) and (2.12) each involve 7 unknowns for 4 moduh and 3 cosines, and are 
not solvable. The missing equation in each group is supplied not by the data but by phase relations 



^LS - ^US - ^LU = {^Lu - - i^Uu - ^5 J - (^L„ - $[/J = (2.14) 

These conditions lead to non-linear relations between the cosines 

COS^i^Ls) + COS^i^US) + COS^($LC/) - 2cos($L5) cos($[/s) COs($Li7) = 1 (2.15) 

COS^(¥l5) + C0S^($[/5) + COS^($L[/) - 2cOs($L5) COs($;75) COs($z,[/) = 1 

Similar relations also hold for the sines. Substituting from (2.11) and (2.12) we get 

al\Sf + ai|L|2 + al\U\'^ - \Sf\Lf\Uf = 2040500 (2.16) 

all's f + al\Lf + al\Uf - = 2040505 

From (2.11) we have 3 equations for moduli 

|5|2 = (oi + 02)-3|L|2 (2.17) 

|C/p = |L|2-l(02 +03) 
|iV|2 = |L|2- 1(02-03) 

Substituting (2.17) into the first equation in (2.16) we get a cubic equation for |Lp = x 

ax^ + bx^ + cx + d = (2.18) 

where o = 3 and 

b = -3[l(oi + 02) + l(a2 + aa)] (2.19) 

c = ^(oi + 02) (02 + 03) + o| + O5 - 3o6 

d = (01+ 02)06 — l(o2 + 03)04 — 2O4O5O6 

From (2.12) we get equations for moduli similar to (2.17), and from the second equation in 
(2.16) we obtain a cubic equation for |Lp = x with coefficients given the expressions (2.18) with 
observables 0^, A; = 1, 6 replacing o^. A; = 1, 6. 

Apart from some factors ^ and the observables aj. and ak,k = 1,6 defined in (2.11) and 

(2.12) are density matrix elements Rep^ + Repy and Rep^ — Repy corresponding to pure initial 
states with target polarization Py = +1 and Py = —1, respectively. The assumptions (2.14) 
mean that the density matrix elements ai..,k = 1,6 and ak,k = 1,6 each form a dccohcrcncc free 
subspace of the reduced density matrix. As the result of a symmetry of the interaction of the 
pion creation process with the environment, these elements decouple from the environment. The 
equations (2.11) and (2.12) will thus not be solved by environment-averaged moduli and cosines 
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TABLE III: Three categories of solutions of Eq. (3.1). 



P < 
V = cos(0) 



F < 
V — cos1l((/)) 



P > 



V = sinh(0) 



ri=+2i?cos(^) ri = i?cosh(|) +iV3i?sinh(|) = i?sinh(| ) + iV3i?cosh(| 



f2 



-2i?cos ^ 



rg = -2i?cos I 



rg = -2i?cosh( f 



= -2i?sinh( I 



but separately by environment-dependent moduli and cosines defined by Kraus representation 45l |. 
On general grounds, there are four interacting degrees of freedom in the interaction of pion cre- 
ation process 7r~p — > 7r~7r"'"n with the environment leading to four sets of moduli and cosines 45, il] 



For any Jmax > 1 the S- and P- wave density matrix elements form a completely autonomous 
subsystem of the reduced density matrix in a sense that it is analytically solvable and determines 
completely the S- and P-wave transversity amplitudes Ar,T = u,d as, well as helicity amplitudes 
An,n = 0,1 with definite naturality [i^. The S- and P- wave density matrix elements of the 
reduced density matrix do not depend on the relative phase uj and involve only the reduced 
transversity amplitudes (2.5). The phase uj is determined analytically in the process of conversion 
of transversity amplitudes At-,t = u,d into helicity amplitudes An with n = 0, 1 nucleon helicity 
non-flip and flip, respectively [3]. As the result the transversity amplitudes (2.6) and (2.7) are 
known up to an overall absolute phase. In the next Section III. we show that there are four 
sets Au{i), Ad{j),i, j = 1,2 of transversity amplitudes solving the system. We shall identify 
these four solutions for S- and P-wave amplitudes with the four environment-dependent S- and 
P-wave amplitudes predicted from the Kraus representation describing the co-evolution of the 
pion creation process with the quantum environment. 



III. DETERMINATION OF REDUCED S- AND P-WAVE TRANSVERSITY AMPLI- 
TUDES FROM MEASUREMENTS OF Rep^ AND Rep°. 

The cubic equations for |Lp and |Lp are just another form of the phase conditions (2.14). The 
equations can be solved analyticaly. To solve for x we write x = y — i/q and require that the cubic 
equation ax^ + bx'^ + cx + d = transforms to the form 

y^ + 3Py + 2Q = (3.1) 

This is accomplished with yo = b/3a and 



Next we define quantities 



There are three categories of of solutions of cubic equation (3.1) [54]. They are given in Table III 



R = sign(Q) vTfl , V = -^>0 (3.3) 
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TABLE IV: Range of momentum transfer —t and dimeson mass m of CERN data on polarized target for 
which high resolution amplitude analyses were performed. 



No. Reaction piat —t m References 













{GeV/c) [GeVlcf 


(JMeV) 




1 


TT 


P - 


TT^TT^ 


n 


17.2 


0.005-0.20 580-1080 


[14, 22] 


2 




P - 


-> TT^TT^ 


n 


17.2 


0.00-1.00 


710- 830 


[12,21J 


3 


7^^ 


'n - 


TT+TT" 


P 


11.85 


0.20-0.30 


360-1040 


^,23] 


4 


7^^ 


- 


-> TT+TT" 


P 


11.85 


0.10-1.00 


710- 830 


^,23] 


5 


TT+n - 


TT+TT" 


P 


5.98 


0.20-0.30 


360-1040 


[16,23] 


6 


7^^ 


- 




P 


5.98 


0.10-1.00 


710- 830 


[16,23] 


7 




- 


K+n 


^P 


5.98 


0.20-0.30 


812- 972 


[18,23] 


8 




- 


K+TT 


^P 


5.98 


0.10-1.00 


842- 942 


[18,23] 



Previous analyses of CERN data on tt p — > vr tt"*"?!, tt"*"?! — > 7r"'"7r p and ii'"'"n — > ii'"'"7r p in 



kinematic regions [17|, Il9|, |2g] which included the study of the complex solutions found in most 
(m,t) bins > and negative values for the third solution X3 = rs — j/Qj irrespective of the 
category of the solution. This allowed later studies 27, 2^, 29] to focus on the category with three 



real solutions. The focus of the present analysis is also on this category. 

We have performed high resolution amplitude analysis of CERN data on vr^p vr^vr+n, 
Tx'^n — > vr+vr^p and K^n K^Tr~p in 8 kinematic regions listed in Table IV. The analyses 
used 1 million of Monte Carlo samplings of error volumes of measured density matrix elements. 
Not all values of density matrix elements correspond to physical amplitudes. The procedure 
produces frequency distribution of physical solutions for the moduli and phases and for various 
observables calculated from the solutions for amplitudes. The average values of the distributions 
for the moduli add up exactly to S = cPa/dtdm and those for the phases satisfy cosine conditions 
(2.15). These averaged values can thus be interpreted as the measured amplitudes. The range of 
distributions determines the asymmetric error bars on the amplitudes and calculated observables. 
Analyses of n^p vr^vr^n using 4 million and 5 million Monte Carlo samplings produce virtually 
identical averaged values and slightly larger errors, indicating the analysis is stable. 

The initial aim of the new study was to resolve resonant strucures in the mass spectra in 
7r~p — > vr^Tr"*"?! at 17.2 GeV/c in the mass range 580-1080 MeV at small momentum transfers 
—i=0. 005-0. 20 (GeV/c)^. The results for the moduli squared of transversity amplitudes in 
7r~p — > 7r~7r"^n are shown in Figures 1-4. We discuss these results and the evidence they present 
for ^"(770) — /o(980) mixing in Section VI.. where we also present additional supporting evidence 
from analyses of several other reactions in Table IV.. 

For each set of solutions for moduli |^(i)|, \A{j)\,i,j = 1,2 we calculate from (2.11) and (2.12) 
the cosines 

\L{i)\\S{i)\ 



COs($Lt/(i)) 



05 



\ums{{)\ 

0-6 

m\\u{i)\ 



13 



Spin Up (+1/2) 



Sp±n Down (—1/2) 



22 
20 
18 

ie 

14 

12 

1 O 



4-3 O 

G 22 

O) 

> 20 

g ie 

I — I 

■ — ■ 14 
1 2 
10 
8 

e 

4 

2 
O 



IL^jl"^ Solution 1 



^H^^- 1 1 h 



H h 



|L 1^ Solution 2 



|L^|^ Solution 1 

.1* 



t t 



t 



|L^r Solution 2 



•I 



0.6 0.7 0.8 0.9 1 1.1 
m ( 71 71 ) 



O.e 0.7 0.8 0.9 1 1.1 

(GeV) 



FIG. 1: Solutions 1 and 2 for P-wave amplitudes and |Ldp for -t=0.005-0.20 (GeV/c)^. 



Spin Up (+1/2) 



Spin Down (—1/2) 



3 

> 2 



O 

e 

5 
4 
3 
2 



Solution 1 



H 1 

|Sul^ Solution 2 



4^ 



ISdl Solution 1 



H 1 1 1 * l *' ^ I 

IS^I^ Solution 2 



0.6 0.7 0.8 0.9 1 1.1 0.6 0.7 0.8 0.9 1 1.1 
m ( 7t:~7t:"^ ) (GeV) 



FIG. 2: Solutions 1 and 2 for S-wave amplitudes jS^I^ and j^dp for -f=0.005-0.20 (GeV/c)^. 
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Ta-irg-et Spin Up (+1/2) Xairget Sp±n Down ( — 1/2) 



e 
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4 






1^ 


Solution 1 






Solution 1 


3 

^ 2 
a> 


















t 








= 1 
^ n 










1 1 




1 1 




1 


It 


1 




j_> <J 

> 








1 1 1 1 
Solution 2 


1 1 


1 1 1 1 
Solution 2 
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FIG. 3: Solutions 1 and 2 for P-wave amplitudes |f/„|2 and \Ud\'^ for -i=0.005-0.20 (GeV/c)^. 

with similar equations for cos{^Ls{j)), cos{^us{j)), cos($l[/(j)). Prom the phase condition 

^L-^S = - ^s) + (^L - ^u) (3.5) 

we obtain 3 equations for sines 

Ca = sin($L5) sin($[/5) = + cos{^lu) - cos{^ls) cos{^us) (3.6) 
Cb = sm{^Ls) sm{^Lu) = +cos{^us) - cos{^ls) cos{^lu) 
Cc = sin($[/5) sm{^Lu) = - cos{^ls) + cos{^us) cos{^lu) 



The system is solvable provided that sign(Cc)=sign(C^CB). Given sin($i:5), the sin($[/5') and 
sin($i;7) can be calculated from Ca and Cb, respectively. With cosine and sines known, the phases 
can be calculated. With 

sin($L5) = eV'l-cos2($£,5), sin($L5) = e^l - cos2($£,5) (3.7) 

where the signs e = ±l,e = ±1. The amplitude analyses assumed positive root for both sines. 
Such procedure is uniqe provided that the phases ^ls and ^ls do not change signs. A change of 
sign would manifest itself as a double zero in the phases $ ls and calculated from the positive 
roots for the sines. In all 8 analyzed reactions in Table IV. we found no evidence of a double zero 
in any of the phases ^ls and $l5 so calculated and apart from the sign ambiguity these phases 
are uniquely determined. 
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Ta-irg-et Spin Up (+1/2) 



Xa.irget Sp±n Down ( — 1/2) 




FIG. 4: Solutions 1 and 2 for P-wave amplitudes \Nu\'^ and |7Vd|^ for -i=0.005-0.20 (GeV/c)^. 

The mass dependence at small —t and the t-dcpcndence at p°(770) mass region of ^ls and 
in Tr~p Tr~TT~^n arc shown in Figures 5 and 6, respectively. Figures 7 and 8 show the mass 
dependence of the phases ^us, ^us sand ^lu, ^lu, repectively. As the figures show, these phases 
are continous and nearly constant functions. Prom this fact we can understand why the phases 
^LS,^LS cannot change sign. Prom (3.6) we see that a change of signs of or ^ls results in 
the change of sign of all phases of the same transversity. This would result in large unphysical 
discontinuities in the phases ^us, ^us or ^lu,^lu- 

In the next Section IV. we show that the phases ^ls and are closely related to the 
positivity of the reduced density matrices p^{Py) = + PyP^ for pure initial states Py = ±1. The 
fact that these phases do not change signs implies that the eigenvalues of the density matrices 
p^{Py = ±1) are all non-zero. 

For negative roots in (3.7) all phases change signs. There are four combinations of the signs 
e and e. The phases with opposite signs correspond to complex conjugate reduced transversity 
amplitudes. For any given solution of moduli |A(i)|, \ A{j)\,i,j = 1, 2 there is a four-fold ambiguity 
in the phases of the reduced transversity amplitudes A = L^U 

L = \L\ exp(ie$L5), ^ = |^| exp(ie$L5) (3.8) 

U = \U\ exp(ie^>i75), U = \U\ exp(ie¥[75) 

In the Section V. we show how the four-fold sign ambiguity can be resolved into a single set of 
phases for each set of moduli |^(i)|, \A{j)\^i,j = 1, 2 in measurements of additional information on 
polarized targets. 
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Target Sp>±n Up> ( + 1/2) 



Tarcget Sp±n Down (—1/2) 
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^ o 
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-r-H 
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PS 



*t*L " *t>s Solution 1 



I I I I |l I I I I 

^t'L " ^i's Solution 2 



*t*L " *t>s Solution 1 



fTTT l|l TTT l|l I T I l|l 1 1 1 

*t*L " *t>s Solution 2 



0.6 0.7 0.8 0.9 1 1.1 0.6 0.7 0.8 0.9 1 1.1 
m(TC^TC*) (GeV) 



FIG. 5: Solutions 1 and 2 for relative phases <i>LS = - and <i>LS = $Ld - for -t=0. 005-0. 20 
(GeV/c)2. 



IV. POSITIVITY AND DECOHERENCE CONSTRAINTS ON REDUCED S- 
WAVE DENSITY MATRIX AND THE PHASES $ls AND 



AND P- 



The measured density matrix elements Re{p^) jj^, and Re{py)'ly are real parts of elements of 



^o^JJ' ^pO^JJ'^^ rjij^g matrices iPu)\xi and iPy)'lx' hermitian and 
so is their measured combination for target polarization component Py 



complex density matrices {Pu!\\ 



yJxx' 



(4.1) 



We will be interested in two specific cases of pure initial states Py = = ±1. Omitting the labels 
^, y, the argument (Py) and suppressing the indices JJ', these reduced matrices have the following 
general form for the S- and P-wave subspace 



/ Pss 

Pis 
POs 



Psl PsO Ps-l \ 
Pii Pio Pl-l 
Poi Poo PO-1 



\ P-ls P-U P-10 P-1-1 / 

Using relations due to hermiticity and P-parity conservation [4 



(4.2) 



\Pk)\'\ 



[Pl) 



Q\JJ' 



-A-A' 



0\J.J'\* 
\\' ) 

X+X'(J)\.JJ' 
\Pk)\\' 



{{Pk)\\ 
i-l 



(4.3) 
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Tax-get Sp±n Up (+1/2) Tajrget Spin Down (—1/2) 




FIG. 6: Solutions 1 and 2 for relative phases ^ls ^ ~ ^s^ and $ls = - <&Sd in P°(770) mass 
region of m=710-830 MeV. 



for k = u,y and similar relations for {p^{Py))\\i , the matrices take the form 

/ Pss Pis POs -P*ls\ 
Pis Pll POl* Pl-1 

Pos Poi Poo —poi 

\ -pis Pl-l -pli Pll ) 



(4.4) 



The positivity of the hermitiam matrices (P requires that det p*^ > 5J, |55|| . Calculations show 
that 



det /9° = (pii + pi-i) {pssPoo{p\\ - P\~\) - '^Ps 



\Po\\ 



(4.5) 



(Pll - Pl-l)\POsf - 2/3oo|/OlsP + 4|/30s||/0ls||/O0l|) 

We split the determinant into two parts det = D + P where 

D = {pii + pi-i){pssPoo{pii - Pi~i) - 2pss{Repoif- 



(4.6) 



{pii - pi-i){RepQsf - 2pQQ{Repisf + ARepQsRepuRep^ 



01 



P= (Pii + /Oi-i)(4|pos||pis||yOoi| -2pss{ImpQiy 



(4.7) 
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Target Spin Up (+1/2) Xairget Sp±n Down ( — 1/2) 
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<I>U " *I'S Solution 1 
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- t|)Q Solution 2 



<|>u - <|>3 Solution 1 



H 1 1 1 |i I I I I 

<t>u - <|>3 Solution 2 
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FIG. 7: Solutions 1 and 2 for relative phases ^us = - *s„ and ^us = ^Ud - ^Sd for -i=0. 005-0. 20 
(GeV/c)2. 

{pii - pi-i){ImpQsf - 2poo{Impisf - ^ReposRepigRepoij 

It is of some interest to note that in general 

det Rep^ = D, det Imp° = (4.8) 

Next we focus on reduced density matrices p^{Py) = P^ + PyP^ with pure initial states = ±1. 
For = +1 we get from the Table I. 

Pss = \S\'', Poo = |i^|', pii - = |?7|2, pn+pi-i = |iV|2 (4.9) 

The expression in the large parenthesis in the term D(Py = +1) is identical to the phase condition 
(2.16) so that D{Py = +1) = as the result of the assumption that the S- and P-wave subsystem 
forms a decoherence free subspace of the reduced density matrix p'j. Using the relations 

Impos = \L\\S\ sin '^>Ls (4.10) 
V2Impis = \U\\S\ sin <^us 
V2Impoi = \L\\U\sin^LU 

and the phase condition (2.15) for the cosines to simplify the term P{Py = +!), the determinant 
for Py = +1 takes the form 

detpO(+l) = \N\^\U\^\L\^\S\\sin^LS? (4.11) 
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Target Spin Up (+1/2) Xairget Sp±n Down ( — 1/2) 




FIG. 8: Solutions 1 and 2 for relative phases ^lu = - ^'i/u and ^lu = ^l^- ^Ud for -t=0. 005-0. 20 
(GeV/c)2. 



We obtain a similar expression for the determinant for Py = —1 

detp°(-l) = \Nf\Uf\Lf\Sf{sm^Lsf 



(4.12) 



The determinants are nonnegative as follows from the positivity of the density matrices p^{Py). 
The density matrices are hermitian matrices and thus can be diagonalized. If all their eigenvalues 
are non-zero, the determinants will be strictly positive. This would imply that the phases 
and cannot change sign, as was actually observed in Section III. from the analysis of the 
experimental data. 



V. RESOLUTION OF SIGN AMBIGUITIES IN THE PHASES OF REDUCED S- AND 
P-WAVE TRANS VERSITY AMPLITUDES IN MEASUREMENTS OF /mp° AND Imp°. 

The S- and P-wave density matrix elements /m/O^ and Imp^ involve interferences between 
reduced unnatural and natural exchange amplitudes of opposite transversity shown in Table II.. 
First we shall show that the elements Im{p^)gl and Im{p^)gl determine the phases un = $jv„ — 
and aN = ^Na~ ^Su of the amplitudes N and N with a two-fold ambiguity for each set of moduli 
\N{i)l \N{j)\,iJ,= 1,2. Prom the Table II. and (2.5) we have 



ri 



V2Im(p°)°}S = Rei-\S\N* + N\S\) = -|5||iV| cosajv + |iV||5| coscktv (5.1) 



si = V2/m(p°)°iS = Im{+\S\N* - N\S\) = -\S\\N\ sinajv - \N\\S\ sinajv 
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where S = (f'a/dmdt. Prom (5.1) we can solve for cosajv and sinajv 

cosaiv = ^^^^(ri + |5||iV|cosajv) (5.2) 
sinajv = ^^^^ (si + |S||]V| sinajv) 



Substituting into cos a^v + sin ajv = 1 we find 



cos ajv , , , — , 
2|5||iV|ri 

where 

|2|c5|2 Icl2il\7i2 / 2 , „2 



^(^-2|5||iV|sisinaAr) (5.3) 



^= |iV|^|5|" - \S\'\NY - {ri + si) 

Substituting into cos^ ajv + sin^a^v = 1 yields a quadratic equation for sinajv with two solutions 

A 



sin ajv 



2\S\\N\{rl + sl) 
where 

4|5|2|iV|2(r2 + s2 



(si±riV5- l) (5.4) 



B 



A^ 

From (5.4) we can now calculate coscFat 
and from (5.3) we obtain 

A 



VS^) (5.5) 



cos ajv 



2|iV||5|(rf + .f) 
sinajv = = — 7, :^(siCibri\/-B — 1) 

where 

A + 2{rl + 4) 



(nC^siV^-l) (5.6) 



C 



A 



Note that for each set of moduli |A(z)|, |A(_7)|, i,j = 1,2 the two solutions for the phases 
'^N{ij),C(N{ij) will depend on the set. 

We shall now show that only one solution is consistent with the remaining equations from the 
Table II. for each set of the moduli and that this solution resolves the four-fold sign ambiguity. To 
this end wc introduce a convenient notation for real and imaginary parts of reduced transversity 
amplitudes A = L,U,N 

A = Ai + iA2 = \A\cosaA + i\A\smaA (5.7) 
A = Ai+iA2 = 1^1 cosSa + «|^| sina^ 
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where the phases and for amphtudes A = L^U arc given by (3.8). The equations for the 
remaining density matrix elements from the Table II. then take the form 

r2 = V2Jm(p°)J}S = LiTVi + LaiVa - LiiVi - L2N2 (5.8) 

52 = V2/m(/9°)J}S = L-iNx - L1N2 + L2N1 - L1N2 

n = -Im{plf\jl = UiNi + U2N2 - UiNi - U2N2 (5.9) 

53 = -/m(p°)i\iS = U2Ni - U1N2 + U2N1 - U1N2 

Suppose that for a given set of moduli |^(^)|, |^(i)| the equations (5.8) and (5.9) are satisfied 
by one of the two solutions for the phases aAr(ij), aAr(ij) for a particular set of signs of phases 
aAi'i),ctAij), A = L,U from (3.8). To be specific, let us suppose this solution for N,N is the 
solution with the positive sign in (5.4) and let us label all amplitudes in this solution A~^,A . 
We now show that the amplitudes N~,N for the other solution with the negative sign in (5.4) 
cannot be a solution of (5.8) and (5.9) for any choice of signs of phases aA{i),olA{j), A = L,U in 
(3.8). 

To prove this statetement let us suppose the contrary and assume that the amplitudes N" , N 
satisfy (5.8) and (5.9) for some amplitudes A", A , A = L,U from (3.8). Since the amplitudes 
A~,A , A = L,U differ from the amplitudes A'^,A'^, A = L,U only in the sign of phases, their 
real parts are the same and the imaginary parts differ at most by sign 

A- = XA+, A2 =XA2, A = L,U (5.10) 

where A = ±1,A = ±1. In the next step we subtract the set of equations with amplitudes N~ ,N 
from the set with amplitudes iV+ , and get a homogeneous set of equations 

LtiN+ - N-) + Lt{N+ -XN^) - L+iNt -Ni) ~ L+(Nt ~ XN2) = (5.11) 
L+(iV+ - AiVf ) - Lt{N+ - N^) + L+(Nt - AiVj - L+(iV+ - iV") = 

Ut (iVi+ -N{) + Ut {N^ - XN2 ) - (Nt - iV7 ) - (N2 - XN2 ) = (5.12) 
C7+(iV+ - AiVf ) - utiN+ - N^) + U+(Nt - AiVn - U+iN^ - ^2") = 

We now write (5.6) in the form 



cos 



^(airiT02Si) (5.13) 
sina^ = |^(aisi ±a2ri) 

and calculate differences and sums 

N+ - ATf = -2a2Si, N+ + = +2airi (5.14) 



A^2^ - = +2a2ri, + TVf = +2aiSi 

essions for amplitudes and N by replac 
in (4.14) with parameters ai,a2 corresponding to the form (4.13) of the equations (4.5) and (4.4) 



We can write similar expressions for amplitudes and N by replacing the paramerters 01,02 
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for coscKiv and sinajv- 



We now examine the equations (5.11) and (5.12) for each possible choice of A and A. The 
case A = A = +1 is excluded as the amplitudes N~,N cannot satisfy the same system as the 
amplitudes N'^,N'^. For the case with A = A = — 1 the equations (5.11) and (5.12) take the form 



si a2 + ai + Lj^a2 — ^2 fli j = (5.15) 

ri ^— L^a2 + L^ai — Lfa2 + -^^2^01^ = 

si(-Uta2 + U2ai + U+a2-U}ai^ = (5.16) 

n (-17^02 + 1^201 - U^a2 + U+ai^ = 

The terms in the parentheses are all different. In particular, the large differences between 
the moduli |L|,|L| and |[/|,|[/| mean large differences in the parentheses for ri at least one of 
which must be non-zero. Since ri has been measured in the CERN experiments on polarized 
targets and is non-zero, the equations (5.15) and (5.16) cannot be satisfied and this case is excluded. 

For the case A = +1 and A = — 1 the equations (5.11) and (5.12) read 

^-Lia2 + L2ai+ Lfa2^si- (^L^a2^ri = (5-17) 
(^-Lia2 + L2ai - Lfa2^ri - (^L^a2^si = 

[-UU2 + + Uta2) si - (u+a2) n = (5.18) 

(-C7+a2 + C7+ai - C/i+a2)ri - (^72+02)51 = 
Combining the first two and the last two equations we obtain two equations of interest 

^-Lia2 + L^ai + 1^0,2^31 = (^-L^ a2 + L2 ai — Lfa2^rl (5.19) 

(-17+ a2 + uUi + t^i+«2) si = (^-UU2 + - U+a2) r\ 

The terms in the parentheses are all different. The parentheses for amplitudes L, L and amplitudes 
[/, V are not proportional to each other since the phases of L, L and C/, \J are 180° out of phase, 
as seen in Fig. 8. Moreover, as we shall discuss in Section VII., the moduli of these amplitudes 
have different widths and structures around /o(980). The equations (5.17) and (5.18) thus 

cannot be satisfied and this case is excluded. The analysis of the case A = — 1 and A = +1 is similar 
with the same conclusion. The solution A^^, selects a unique set of phases of amplitudes L, \J 
and L, C/. The equations (5.8) and (5.9) change when the phases of these amplitudes change sign. 
The change results in a different solution for amplitudes N, N which is not compatible with the 
data in (5.1). 

We conclude that the measurements of iS-and P-wave density matrix elements Imp^ and /m/3° 
listed in Table II. unambiguously select a unique solution for S- and P-wave reduced transversity 
amplitudes. In effect, for any set i,j of solutions for the moduli there exists only one set of signs of 
phases in (3.8) such that the solution of equations (5.8) and (5.9) for Ni, N2, Ni, N2 satisfies the 
conditions |A''(z)p = A^^ + and |iV(j)p = + A^2- These measurements were not feasible in 
the 1970's when the first CERN measurements on polarized targets were done. With the advent 
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of advanced frozen spin targets the direction of target polarization can be selected at will with 
high degree of polarization [55|. The next generation of measurements of pion creation processes 
on polarized targets will therefore provide data on the complete reduced density matrix. 

VI. EVIDENCE FOR p°{770) - /o(980) MIXING IN w-p TT-w+n AND n+n n+Tr-p. 

We first recah the equations (2.17) for the moduh jS'p, and \N\'^ 

\S\^ = (ai + a2) - 3\L\^ (6.1) 
|C/|2 = |L|2-l(a2 + a3) 

\N\' = \Lf - \{a2 - a,) 

The data components ai + a2, 5(02 + 03) and 5(^2 — 03) were calculated from those Monte Carlo 
samplings of density matrix elements in the error volume of the data for which physical solutions 
for the amplitudes were found. The results for ai + a2 and ^{0,2 + ^3) are shown in Figure 9. The 
results for ^(02 — 03) are very similar to those of ^(02 + 03) and are not shown. The Figure 9 
shows pronounced /?°(770) peaks for target spin "down" components and a suppression of ^"^(770) 
in the target spin "up" components. There is a pronounced dip at 970 MeV in target spin "down" 
components and a dip at 1030 MeV in target spin "up" components corresponding to the /o(980) 
resonance. 
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FIG. 9: Data components ai + 02 of the moduli |S'tP and (02 + as)/^ for the moduli |?7rP, t ^ u,d for 
-i=0.005-0.20 (GeV/c)2. 
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Figure 1 shows the two solutions for the P-wave amphtudes and l^dp- Remarkably, 

the two solutions of the two cubic equations reproduce closely the general features of the data 
components in the Figure 9. Figure 2 shows the two solutions for the 5- wave amplitudes \Su\'^ 
and IS'dp. The rho-like resonance and /o(980) are clearly resolved in both solutions for jS'dP, 
in contrast with our previous low resolution amplitude analyses based on 40,000 Monte Carlo 



samplings [27|, |28|, l29|]. This peak at 770 MeV with a width at half-height of 155 MeV has its 
origin in the p^{770) peak of the data component ai + 02 which survives the subtraction of 
3|L(ip. Just like the p^{770) peak in data components ^(02 + 03) and ^(02 — as) for target spin 
"down" contributes to p^{770) peaks in the amplitudes \Ud\'^ and jA^'^P (see Figures 3 and 4), 
so the p°(770) peak in the data component ai + 02 gives rise to the presence of p^{770) in the 
amplitudes 15'^^. Similarly we can see that the suppressed p^{770) peak in the data components for 
target spin "up" propagates to all target spin "up" moduli | ^wp, including 5- wave amplitude 

The data components 04, 05 and ag provide independent evidence for the presence of p^{770) 
in the 5- wave amplitudes. Figures 5 shows relative phases ^ls = ^L„ — ^s„) ^LS = ~ ^Sd- 
For target spin "up" the phase ^ls is constant and small, reflecting the similarity of \Su\'^ and 
For target spin "down" the phase ^ls is near zero in the Solution 1 and small and slowly 
varying in Solution 2 below 1000 MeV. These results indicate that the amplitudes Sr and Lr, 
T = u,d are in phase. Since the amplitudes Lt resonate at p^{770), so must the amplitudes 5,-. 

This conclusion is in agreement with the phases ^us = ^Uu ~ ^Sui ^us = ~ ^Sd shown 
in Figure 7. These phases are nearly constant below 1000 MeV and indicate that the amplitudes 
Sr and Ur are 180° out of phase. Again, since Ur resonate at ^"^(770), so must the amplitudes Sr- 
The Figure 8 shows that the relative phases of the amplitudes Lr and Ur are nearly constant and 
180° out of phase, as expected from the two resonating amplitudes. 

The observation of p°(770) - /o(980) mixing in the 5-wave amplitudes brings up the question 
of p^{770) — /o(980) mixing in the P-wave amplitudes. Figure 1 shows a clear dip at 970 MeV 
in the amplitude l-L^^p followed by a sharp rise. The dip occurs at the mass of scalar resonance 
/o(980) and is more pronounced in the Solution 2. Figure 3 shows unexpected resonant structure 
(a bump) at the /o(980) mass in the Solution 1 of the amplitude \Ud\'^- These findings are a 
strong indication for p^{770) — /o(980) mixing in the P-wave amplitudes. Since /o(980) resonance 
is produced by unnatural exchange, we expect /o(980) to be suppressed in natural exchange 
amplitudes \Nr\'^ as evidenced by the data in Figure 4. 

Independent evidence for yo'^(770) — /o(980) mixing in the S'-wave amplitudes comes from the 
amplitude analysis of the CERN data on the t-dependence of the density matrix elements in 
the ^^^(770) mass region (reaction 2 in the Table IV.) for -t=0. 0-1.0 GeV/c^. In Fi gure 10 we 
compare the results for the amplitudes l/SrP and l^rP- We find that the moduli have similar 
shapes in all solutions for both target spin "up" and "down" amplitudes. This will happen if the 
two amplitudes both resonate at p^{770). 



This conclusion is supported by the results for the relative phases between amplitudes Sr 
and Lr in Figure 6 and amplitudes Sr and Ur in Figure 11. The phases ^ls = ~ ^Su ^^^d 
^LS = ^Ld ~ ^Sd small and nearly constant below -v/— i ^ 0.6 GeV/c, indicating that the 
amplitudes Sr and Lr are nearly in phase. The phases ^us = ^Uu ~ ^Su ^us = ^Ua ~ ^Sd 
are also nearly constant and the amplitudes Sr and Ur are nearly 180° out of phase. Since the 
amplitudes Lr and Ur resonate at p°(770), so must the amplitudes Sr- 
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FIG. 10: Solutions 1 and 2 for the moduli |S'tP (solid lines) and \Lr\'^ (dashed lines), t ^ u,d in the p'~'{770) 
mass region m=710-830 MeV. 

The analyses of CERN data on Tr~^n yr+vr^p on polarized target at 11.85 and 5.98 GeV/c 
confirm the expectation of p^{770) — /o(980) mixing at other energies (reactions 3 and 5 in Table 
IV.). Figures 12 shows the moduli \Su\'^ and {Sdl"^ of the S-wave amplitudes at 11.85 at momentum 
transfers 0.20 < —t < 0.40 (GeV/c)^. Despite low statistics, the mass spectra show evidence for 
p^{770) in the 5- wave amplitudes also at these energies. 

There is no P-wave and isospin 1 = 1 component in the tt^tt^ states in Tr~p — > vr'^vr'^n process. 
We thus do not expect the p^{770) peak to appear in the S'-wave intensity in this reaction. The 
p^{770) interpretation of the observed rho-like resonance in the S-wave in ■K~p 7r~7r~^n thus 
explains the lack of evidence for a narrow rho-like scalar resonance in tt~p — > vr'^vr'^n first observed 
at CERN in 1972 [3] and recently confirmed in high statistics measurements at 18.3 GeV/c at 
BNL 0. 



VII. TEST OF ROTATIONAL AND LORENTZ SYMMETRY IN Tr p ^ TT Tr+n. 

It is generally expected that the position and the width of the p^{770) peak as observed in 
the spin averaged cross-section will be faithfully reproduced on the level of spin amplitudes. The 
CERN data on polarized target show that this is not the case. 

From Figures 1,3 and 4 we see that the p^ production is suppressed in all target spin "up" 
amplitudes while the target spin "down" spectra dominate the p^ production. The width at 
half-height of the peaks in the longitudinal spectra \Lu\'^ and l-L^p is the expected 155 MeV. 
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FIG. 11: Solutions 1 and 2 for relative phases $£/s = — $s„ and $(75 = ~ in the p°(770) mass 
region m=710-830 MeV. 

However, the width shows dilferent values in different transverse spectra. In the spin "down" 
spectra the width is narrower at 120 MeV in \Ud\^ but wider at 180 MeV in |iVdP- These values 
are reversed in the spin "up" spectra with 180 MeV and 120 MeV in amplitudes jf/^p and lA^^^P; 
respectively. Such large variations in the width are entirely unexpected and appear anomalous. 

The rotational symmetry of strong interactions requires that the width and mass of a resonance 
do not depend on its helicity. It also prevents mixing of scalar and vector resonances in the same 
partial wave with a definite spin J. But what if rotational symmetry is broken in pion production? 
Then the widths of p^ in the P-wave amplitudes L^, Ud and could differ and p^{770) — /o(980) 
mixing could occur in both S- and P-wave amplitudes. 

But how realistic is the hypothesis that rotational invariance is violated in pion production? 
The breaking of rotational symmetry implies a breaking of Lorentz symmetry. The breaking 
of Lorentz symmetry has been recently examined theoretically and experimentally in several 
different contexts. In 2001, Gale and collaborators at McGill examined vector - scalar mixing of 
/9'^(770) — ao(980) resonances in dilepton production in hot and dense matter. The medium-induced 
breaking of rotational and Lorentz symmetry leads to mixing of different spin states even when the 
interaction Lagrangian respects all required symmetry properties [31]. While the effects of broken 
Lorentz symmetry are large at high temperature T, they vanish for T — > 0. Their work inspired the 
idea to test rotational and Lorentz symmetry in CERN data on ir^p ir^ir'^n on polarized targets. 

The transverse amplitudes Ur and Nj- which show very different p^ widths are a mix of dipion 
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FIG. 12: Solutions 1 and 2 for the moduli \Sr\'^, T = u,dm n+n Tr+n-p at 11.85 GeV/c and -4=0.20-0.40 
(GeV/c)2. 



+1 and -1 helicities and thus are not suitable to test the rotational symmetry. To test the rotational 
symmetry we need transversity amplitudes with a definite dipion helicity . The required amplitudes 
are H^^,H~^^ Hj^^,H^^ and they are related to the amplitudes Ur and 



(7.1) 



Their partial wave intensities and polarizations can be calculated from the data on polarized target 



(7.2) 



P(i?„) = -2Re{H+'H-'*) = - |Ar„|2 

The equations for I{Hii) and P{Hd) are similar. For with zero helicity we have longitudinal 
amplitudes and L^. It is convenient to relabel them as 



Their moduli squared shown in Figure 1 are then the longitudinal intensities 



(7.3) 
(7.4) 
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FIG. 13: Solution 1 for the intensity I{Ht) and the interference P{Hr), r 
amplitudes H~^^ and with definite dipion heliciti A = ±1. 



u,d of the P-wave transversity 



Figures 13 and 14 show the transverse intensities I{Hu),I{H(i) and the polarizations 
P{Hu), P{Hd) for the two solutions. The intensities show a clear single peak with the same width 
at half-height of ~ 155 MeV for both target spins in both solutions. This transverse width is 
exactly the same as the longitudinal width in the intensities I{H^) and I{H^). This indicates 
that the p° poles in all helicity amplitudes H^,H^, A = 0, +1, -1 have the same width ~ 155 
MeV. 

This conclusion is supported by the observation that the polarizations P{Hu) and P(Hd) 
are both small and vanish at the peak mass in I{Hu) and I{H(i), respectively. The pairs of 
amplitudes {H^^,H~^) and (H^^ , H^^) with opposite transverse helicities are 90° out of phase 
at the peak mass and move around 90° at other masses. Such behaviour is consistent with the 
two amplitudes in each pair being dominated by the same /9° Breit-Wigner pole interfering with a 
weak non-resonating background. 

On this basis we can conclude that there is no evidence for any dependence of p^ width on 
its helicity from the CERN data on pion production on polarized target. To explain the large 
differences in the p° witdth observed in the spectra and jA^dP we note from (7.2) 



\Ud\' = -iim+pm) 



(7.5) 
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FIG. 14: Solution 2 for the intensity I{Ht) and the interference P{Hr), t = u,doi the P-wave transversity 
amplitudes i?"*"^ and with definite dipion heliciti A = ±1. 

Figures 13 and 14 show that the polarization P{Hj) is negative and broad in the mass region. 
As the result of the opposite signs of P{Hd) in (7.5), the p° width is narrower in |J7dp and wider 
in \Nd\^. This apparent difference is entirely due to the interference of the amplitudes H'^^ and 
H'^^. This interference arises from the presence of a non-resonating background (continuum) in 
these amplitudes. The relations for the target spin "up" spectra \Uu\'^ and \Nu\'^ have the same 
form as Eqs. (7.5). However, now the polarization P{Hu) is positive and broad around the 
mass and it has therefore an opposite effect on their widths, as observed. 

Rotational invariance of strong interactions and the particle interpretation of resonances also 
require that the p^ mass rup does not depend on its helicity A. The test of this aspect of rotational 
symmetry in pion production is somewhat less direct since the positions of the observed peaks in 
the P-wave spectra need not correspond exactly to the p^ mass as the result of interference effects 
with non-resonating background. 

Figures 13 and 14 show that the transverse intensities /(i?d) and I{Hy) for target spins "down" 
and "up" peak at 770 and 790 MeV, respectively. In Figures 1 we see that the longitudinal 
intensities I{H^) and I{H^) for target spins "down" and "up" peak at 790 and 810 MeV, 
respectively. The apparent dependence of the peak positions on the target spin can be understood 
as the intereference of the p^ pole with the target spin dependent background. We may conclude 
that there is no evidence for the dependence of the p^ mass on its helicity from the CERN data 
on pion production on polarized target. 




30 



The amplitude analysis of tt~p vr^-Tr+n on polarized target presented in Section III. provides 
the first direct experimental evidence for the independence of the mass and width of a resonance 
on its helicity expected theoretically from the rotational symmetry of strong interactions. If strong 
interactions do not violate Lorentz symmetry in pion production then we have to look somewhere 
else to find an explanation for the observed p°(770) — /o(980) mixing. 



VIII. NON-EXISTENCE OF A UNIQUE ANALYTICAL SOLUTION AND 
THE CENTRAL HYPOTHESIS. 



Apart from a few Monte Carlo samplings in one m-bin in reaction 1 and in 4 f-bins in reaction 

2, all Monte Carlo samplings of the data error volume that yield physical solutions for the moduli 
and cosines require that R > 0. This numerical result means that the two solutions for the moduli 
have analytical form 

\Lr{l)f = -yo + Rcos{^) + VSRsm{^) (8.1) 

|L,(2)|2 = -yo + i?cos(^)-A/3i?sin(^) 

where we suppressed the subscript r on r.h.s. of (8.1). It follows from (2.17) and (6.1) that the 
two solutions for the moduli \St-\'^, \Ur\'^ and \Nt\'^ have a similar analytical form. 

Prom the Table III. we see that the condition + = defines a surface Eo,t which divides 
the error volumes of observables a^, and a^, k = 1,6 into two separate regions and S-,r 

with + > and + P^ < 0, respectively. In S+,t there are two complex conjugate 
solutions and a negative real solution for [L^-p. In E_^t- there are two distinct positive solutions 
and one negative solution for |LrP- On the boundary So,r there are two degenerate solutions 
corresponding to sin(^) = and comprising one positive solution, and one negative solution. The 
two conditions + P^ = for the observables ak, and a^, k = 1,6 divide the data error volume 
of the measured density matrix elements similarly into two regions S+ and S_ characterized by 
complex conjugate and positive solutions for the moduli, respectively. 

From the standard Quantum Field Theory we expect a unique physical solution. This 
expectation raises a question whether such unique physical solution could be identified with the 
single positive solution on the boundary Sq. This solution would describe a unitary evolution 
from pure initial states to pure final states, would be time-reversible and would represent the pion 
creation process as an isolated event in the Universe. 

We investigated this possibility numerically as follows. The difference between the solutions, 
which we call level splitting, 

Sr = |ir(l)|^ - |ir(2)|^ = 2 Rr smfy (8.2) 

o 

measures the distance of the solutions from the boundary Sq.t- We looked at a number of physical 
solutions constrained progressively closer to the boundary by imposing a condition on the solutions 
that 6t < S where the upper bound S was progressively reduced. The amplitude analysis was done 
with normalized amplitudes such that 

'^'^ E \Au\' + \Aa\' = l (8.3) 



dmdt ^ c, r TT,., 

A=S,L,U,N 
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TABLE V: Total number of physical solutions in progressively closer vicinities of the boundary Eq. The 
parameter S is upper limit on level splitting for normalized amplitudes \Lr\'^,T — u,d. For S — free results 
were scaled from sample size 5.10^ to 100.10^. 



6 free 


0.0010 


0.0005 


0.00025 


0.0001 


Sample size 20(5.10*^) 


100.10^ 


100.10^ 


100.10^ 


100.10^ 


Physical solutions 66008240 


12499 


1075 


82 


2 



For each limit 6 the analysis was done with 100 million Monte Carlo samplings of the data error 
volume. The results are shown in Table V. 

The results clearly show that the number of physical solutions dramatically decreases as the 
boundary is approached. The real solutions on the boubary are thus not physical. Physical 
solutions are confined to a region away from the boundary leading to two distinct solutions for 
the moduh \A{i)\'^ = \Au{i)\'^ and = \Ad{j)\'^, i,j = 1,2. This multiplicity of solutins is 

inherent in the data in all 8 analyzed reactions listed in the Table IV. 

We have shown in Section V. that the complete measurements on polarized target uniquely 
determine the phases of reduced S- and P-wave transversity amplitudes A and A, A = S, L, U, N. 
In a sequel paper we shall show that the relative phase uj in the S- and P-wave transversity 
amplitudes A^- is uniquely determined by a process of conversion of transversity amplitudes 
into the helicity amplitudes. Up to an inessential overall phase, the S- and P-wave transversity 
amplitudes Au{i) and Ad{j) are thus uniquely determined for each of the four solutions i,j = 1, 2. 

The four solutions cannot be distinguished in measuremnts of reduced density matrix on 
polarized targets since for all measured density matrix elements = p^f{ij) for all However 
the predictions for full density matrix, and thus for the measurable recoil nucleon polarizations, 
are all different pf{ij). We could hope that the measurements of recoil nucleon polarization would 
select one valid physical solution. But we are also free to think that all four solutions are valid 
physical solutions and that their multitude reflects new and important aspects of reality. 

This is our central hypothesis. It implies that the measured final state density matrix is a mixed 
state of the four physical solutions 

Pf{e(t>, p) = piip/(ii, e(t>, p) + pi2P/(i2, e(t>, P) + p2ip/(2i, ^0, p) + ^22^/(22, ^0, p) (8.4) 

where the probabilities pij satisfy condition ^ pij = 1 and P is the target polarization. In 

a sequel paper [46(] we show that each set Au{i), Ad{j),i, j = 1,2 of transversity amplitudes 
uniquely determines the 5- and P-wave subsystem of the density matrix pf{ij,6(p,P). The 
transversity amplitudes for higher spins J also come in four sets A-l^{i),Al^{j),i,j = 1,2 so that 

the corresponding density matrix pf(ij,9(p,P) is fully determined. The probabilities pij can be 
experimentally determined in measurements of recoil nucleon or hyperon polarization on polarized 
target [ij ]. 

The hypothesis (8.4) implies that pure initial states evolve into mixed final states in the 



pion creation process 7r~p — > 7r~7r~''n. In our previous paper 45|] we have shown how the 
conditions for unitary evolution of pure initial states into pure final states are violated by data 
in TT~p 7r~7r"*"n, implying a non-unitary evolution in pion creation process. In quantum 
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theory such 
environment 



non-unitar\ 

4^, H m 



evolution arises in interactions of open quantum systems with an 
Physically the hypothesis (8.4) thus means that the pion creation 



process can be thought of as open quantum system interacting with a quantum environment. 
Such interaction is described by a unitary co-evolution operator U replacing the standard 5-matrix. 



In the next Section IX. we show how the analytical solutions for the transversity amplitudes 
reveal the existence of a new quantum number g characterizing the quantum states of the envi- 
ronment. In Section X. we show that these analytical solutions can be identified with co-evolution 
amplitudes which are matrix elements of the co-evolution operator U corresponding to the quantum 
states of the environment. 



IX. LEVEL SPLITTING AND THE QUANTUM NUMBER g. 

When we look at the mass spectra of the amplitudes |C/p and \N\'^ in Figures 1-4 we 

notice that the two solutions are very close. The spectra carry definite quantum numbers of the 
target nucleon transversity corresponding to a two-level quantum system. They also carry the 
angular quantum numbers J, A of the two-pion state corresponding to one-level and three-level 
quantum systems for the S- and P-wave spectra, respectively. If we think of the mass spectra 
as variable spectral lines in dipion mass, the two solutions suggest a level splitting due to some 
unknown interaction associated with a new two- valued quantum number g = ±1 corresponding to 
the signs ±1 in the analytical solutions (8.1). From the hypothesis (8.4) it follows that the level 
splitting arises from the interaction of the pion creation process with a quantum environment. 

The difference between the two solutions for the P-wave amplitudes has the same value 

A, = |L,(1)|2-|L,(2)|2 (9.1) 

indicating that the level splitting does not depend on the dipion helicity A. It depends on the spin 
J and for the S'-wave 

\Sr{l)\^ - \Sr{2)\^ = -3Ar (9.2) 

The results for A,- are shown in Figure 15. The largest level splitting occurs in A^ for masses in 
the vicinity of the /o(980) resonance. 

We see from (8.1) that the P-wave moduli have a general form 

\Ar{l)\^ = {XArf + {Yrf (9.3) 
\Ari2)\^ = {XArf-(Yrf 

where (i^)^ = ^A,-. Next we define amplitudes associated with the quantum number g-j- 

ZaAQt) =exp{igrXAr)\Aril)\ (9.4) 
where the phase Xat is given by the level splitting 

tan A^^ = — ^ (9.5) 

-X-At 

The intensities (9.3) then read as interferences 

\Ar{l)f = Re{ZAr{+l)Z*A,i+l)) (9.6) 
|A,(2)|2 = Re{ZAri+l)Z*AA-i)) 
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FIG. 15: Level splitting A,-, i.e. the difference between the Solution 1 and Solution 2 for the unnormalized 
moduli squared |L^|^, r = m, d. 

In the first interference term the pliase difference {+1)Xat — (+l)AAr = 0. Tlie phase difference in 
the second interference term is non-zero {+1)Xat — {—^)Xat = 2(-|-1)A^t- so that 

|A(2)|2 = cos(2AAr)|A,(l)|2 (9.7) 

Next we associate the solutions for the P-wave amplitudes with the quantum number 

Aril) = Ari+1) = 1^^(1)1 expi$Ar(l) = exp(+iAAr)|^r(l)| exp(i^'Ar(l) (9-8) 

Ar{2) = Ar{-1) = |^r(2)|expi$Ar(2) = exp(-iAAr)|^r(2)|exp(i*^^(2)) 
With '^At{+^) = *At(1) and '^At{—^) = *Ar(2) we can write the amphtudes ^r(5) in the form 

^,(+1) = ZaA+I) cxp(i^Ar(+l)) (9.9) 



Ar{-1) = VcOs(2AAr)^Ar(-l) exp(z*Ar (-1)) 



= VcOs(2AAr) exp{-2iXAr)ZAr{+'i-) exp(?*Ar(-l)) 

With these expressions we find 

lAri+l)^ = Re{Ari+l)AU+l)) = Re{ZAr{+l)Z*A,i+l)) 



(9.10) 
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= Re{Ar{-l)Al{-l)) = Re{ZAr{+l)Z*Ar{-i)) 
For the S-wave we have [^^(l)^ < |S't-(2)P so we have to define 

ZsAgr) = eMi9TXsr)\SA2)\ (9.11) 

Then the moduli are interferences 

|5.(1)|2 = Re{ZsA-l)Z*sA+^)) (9.12) 
|S,(2)|2 = Re{Zsr{-l)Z*sA--^)) 

and the amphtudes associated with quantum number read 

Sr{l) = Sr{+1) = VcOs(2A5r)^5r(+l) exp(i*5r (+1)) (9.13) 
5^(2) = Sr{-1) = Zsr{-l)eMi^Sr{-l)) 

where *5r(+l) = *5r(l) and ^'5r(-l) = *Sr(2). 




FIG. 16: Level splitting phases Xat for amplitudes A = S,L, U, N. Maximum value of the phases is 45°. 

The level splitting phases Xat are shown in Figure 16. From the definituion (9.5) we see that the 
phase Xat describes the relative strength of the interaction with the environment in the amplitude 
At-. For small Xat the interaction has a small relative effect. The effect is maximum at 45". when 
|^T-(2)|^ = 0. The interaction has a larger relative effect in the 5- wave amplitudes Sr than in the 
P-wave amplitudes Lr with zero dipion helicity. The effect is about the same in the transverse 
amplitudes U-j- and Nt- The largest overall effect is in the A^ amplitudes near the /o(980) resonance 
mass. A characteristic feature of all level splitting phases are phase fluctuations. 
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X. QUANTUM STATES OF THE ENVIRONMENT. 



In order to introduce the central concept of co-evolution amplitudes we first briefly review the 
Kraus representation for for reduced density matrices of open quantum systems interacting with 
an environment. It is the co-evolution amplitudes which involve the interacting degrees of freedom 
of the environment and with which we shall identify the solutions for the transversity amplitudes. 

The co-evolution of an open quantum system S with a quantum environment E' is a unitary 
evolution [48] 

Pf{S, E) = UpiiS, E)U+ = Up^{S) ® p^{E)U+ (10.1) 

The initial state of the environment is in general a mixed state 

Pi{E) = ^pee'\ee >< ee\ (10.2) 
t 

where |e£ > are quantum states of interacting degrees of the environment and = 1. The 

Hilbert space of the environment has a finite dimension. It is given by a condition (\.\m.H{E) < 
dimHi{S) diniHf{S) [4^. After the interaction the system S is fully described by reduced density 
matrix given by Kraus representation 

Pf{S) = TrE{pf{S,E)) =Y,T.P-^-^^rnP^{S)S+ (10.3) 

£ m,n 

where the operators S^m =< Gi\U\em > satisfy a completness relation Yl S^eSim = I- 

I m,n 

In our next step we associate the two solutions for transversity amplitudes Au{i) and A(i{j), 
i,j = 1,2 with two single qubit states \i > and \j >, respectively. Then the hypohesis (8.4) 
allows us to identify the four degrees of freedom of the environment | > allowed by the condition 
dim.H{E) < dimHi{S) dimHf{S) = {2sp + l)(2s„ + 1) = 4 with the four two-qubit states \ei >= 
\i > \j >. Since the transversity amplitudes can possess only one solution at a time, the co-evolution 
amplitudes A^'J_(£m) 

Aj;^^{£m) =< JA7?,r„| < ee\U\em > \OTp > (10.4) 
must be diagonal for any dipion spin J and naturality r] 

Ailiim) = Ai'^^iii)6,m = Aj^^^iij, ij)5ij^,j, ^ 4' (^J^^.'.^'i' (10-5) 

where 

^ili^j) =< JXv,rn\ < ij\U\ij > \0u >= Ailii) (10.6) 

In (10.4)-(10.6) the state |Orp > is the Tr~p initial state with pion spin and proton transversity 
Tp = u,d. In (10.6) A'!j'^{i), Aj^^^{j),i,j = 1,2 are solutions for transversity amplitudes with dipion 
spin J and helicity A generalizing the S- and P-wave amplitudes Au{i),Ad{j)- 



36 



Instead of using the solution qubits \i > and | j > to define the quantum states \i > \j > of the 
environment and the co-evolution amplitudes A'!j'^{ij), we could have used the qubits \gu > and 
\gd > to define equivalent states \gu > \gd > of the environment and the equivalent co-evolution 
amplitudes A'^^^{gugd) = A'^^^{gr). These states more closely reflect the qubit nature of the 
interacting degrees of the environment and may possess a deeper physical meaning. Recall that 
the solutions for the moduli are distinct and well away from the boundary Sq. This suggests that 
the quantum numbers gr are good quantum numbers suitable to describe the interacting degrees 
of freedom of the environment. 



XI. NON-UNITARY DYNAMICS OF /(770) - /o(980) MIXING. 



Non-unitary evolution to mixed final states (8.4) in 7r~p yr^vr^n arises from a CPT violating 
interaction of the pion creation process with a quantum environment. The hypothesis of the 
existence of such an environment explains naturally the observed p''(770) — /o(980) mixing. We 
can see how this may work using the following mechanical analogy. Consider a pendulum of mass 
m and length L oscillating with a natural frequency w, and an ensemble of several other penduli 
with various masses rrii and lengths Lj. When the other penduli are isolated from the oscillating 
pendulum they do not oscillate. However, when all the penduli are attached to a common rod all 
penduli begin to oscillate with the same resonant frequency although with different amplitudes. 
The interaction of the penduli with a common environment - the rod - allows a resonance from 
one pendulum to "leak" into the other penduli which have different natural frequencies. Similarly 
we can imagine that the resonant qq modes produced in ■K~p TT~TT~^n process interact with an 
environment which allows resonances to "leak" into different two-pion partial wave amplitudes. 



The requirement that the Kraus representation leaves invariant the spin formalism used in the 
data analysis necessitates j45^ that the co-evolution amplitudes 

^xlriar) =< JXv, rn\ < gu9d\U\gugd > |0t > (11.1) 

transform under P-parity as a two-body P-parity conserving process vr^ +p " J(7r^7r^)" +n with 
parity P = (—1)'^ for the dipion states " J(7r~7r+)". This means that there is no vector associated 
with the quantum states \gugd > of the environment. In particular, there is no energy-momentum 
associated with these quantum states. The interacting hadrons conserve their energy-momentum 
and there is no exchange of energy-momentum with the environment, in agreement with the original 
proposal by Hawking for particle processes interacting with quantum fluctuations of the space-time 
metric 33|. Instead, the interaction with the environment is a non-dissipative dephasing process. 



The co-evolution amplitudes U fi can then be written in a form 

Ufi = If^ + i{27r)^6\Pf - Pi)Tfi (11.2) 

where Pi and Pf are total four-momenta of the initial and final hadron states and Tfi is the 
transition matrix for the process \Tr~p > +\gugd >^ Ivr^vr+n > +\gugd >■ 



To construct a model of interaction of the pion creation process S with the environment E we 
can think of the transition operator T as a product of three evolution operators 

T = TfETpTi (11.3) 

corresponding to three different stages of the pion creation process. In the first stage the operator 
Ti maps the initial hadron state \i >= |0r > from Hilbert space Hi{'ir^p) to a vector \iint > in 
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an intermediate (hidden) Hilbert space Hint[qqn) of resonant and non-resonant (gg) modes and 
neutron states wliieh form an orthonormal basis in Hint{qqn). In general, these states could include 
formations also of four-quark states such as diquark-antidiquaurk states. For simplicity we assume 
that the resonant modes correspond to quark-antiquark states in LS'-coupling. Focusing on the 
resonant modes component of \iint > we have 

Ti\i > \gu9d >= \iint > \gu9d >= (11-4) 
< K ix{qq)7], Ik ,TK^i\Ti\QT > \Kn{qq), Ik,tkh >)|5'«5'd > 

where K and /x are the spin and helicity of the qq state, Ik is its isospin, tk^ is the transversity 
of the recoil neutron in the intermediate process 7r~p — ^ qqn and r/ is its t-channel naturality. The 

resonant states propagate with Breit-Wigner amplitudes a Kirn?) and the coherent state \iint > is 
modified in Hint{qqn) by the action of the propagation operator Tp 

\iint,p > \gu9d >= {Tp\iint >)\gugd >= (ll-5) 

< Kn{qq)r],lK,TKn\Ti\QT > Tp\KiJ,{qq),lK,TKij. >) bwfd >= 

< Kfj,{qq)ri,lK,TK,,\Ti\OT > aK{m?)\Kfj,{qq),lK,TKij. >) > 

As the resonant modes propagate they interact with the environment. It is crucial to realize 
that the interaction involves simultaneosly the whole coherent state \iint,p >■ In this final stage 
the the states [Kn^qq), lK,TKf_i > \gugd > are mapped into the Hilbert space IIf{n~7r'^n) (gi H{E) 

TfE\Kn{qq),lK,TKn > >= (11-6) 

X] < 9ugd\ < JK'^~'^^),lK,Tn\TfE\Kn{qq),lK,TK^ > \gugd > |JA(7r~7r+),/i^,rn > \gugd > 
JX 

where we have assumed isospin conservation in transitions from \K fi{qq), Ik tTk^i > to 
I JA(7r~7r+), /x, T„ >. The total co-evolution amplitude will involve isospin amplitudes 

< JA(7r-7r+)r/, Ik, Tn\ < gugd\T\gu9d > |0r >= (11.7) 
^ < JA(7r"7r+)77,/K,rn| < gu9d\TfE\9ugd > \Kn{qq),lK,TKix > 

aK{m?) < Kji{qq),lK,TKij\Ti\QT > 

where Ik = /(tt^tt"*") = 0, 1 is the isospin of the it~it~^ state. The non-resonant modes with Ik = 
give rise to a non-resonating amplitude < JA(7r^7r+), /(7r^7r+) = 2,rn| < gu9d\T\gugd > \^'t > 
given by an expression similar to (11.7) with a replacements ai^(m^) — > 1 and gg qqqq- 
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In standard spectroscopy there is no environment interaction and the matrix elements of TfE 
are diagonal with J(7r~7r"'") = K{qq) and A(7r~7r+) = fJ-{qq). As a result, there is no mixing of reso- 
nances in two-pion partial wave amplitudes as the basis vectors in Hint{qqn) are mapped into basis 
vectors in Hf{'iT^7r^n). The effect of the interaction of the resonant modes with the environment 
are transitions from the resonant modes \Kn{qq, Ik,tk^ > to two-pion states \JX{tt~it~^), lK,Tn > 
as the basis states in Hint{qqn-) are mapped into superpositions of basis vectors in Hf{TT^'iT~^n). 
The non-diagonal matrix matrix elements of TfE allow for the resonances to "leak" into different 
two-pion partial wave amplitudes and account for the observed p^{770) — /o(980) mixing. The 
interaction with the environment is non-local and CPT violating which accounts for the violation 
of the conservation of angular momentum and parity in the non-diagonal transitions between these 
two different quantum systems \^ . 



XII. ENTANGLEMENT OF tt 7r+ ISOSPIN STATES AND CPT VIOLATION. 



Before we can write down the final form of the co-evolution amplitudes we need to clarify 
the isospin structure of the produced two-pion states. The Generalized Bose-Einstein symmetry 
assumes that all particles in an isospin multiplet are identical particles. This symmetry implies 



that spin and isospin of two-pion non- interacting states must satisfy condition J + I = even [53|] . 
As the result, the symmetrized 7r~7r+ isospin state is symmetric for J = even and antisymmetric 
for J = odd 

J=even |5 > = ^(Ivr" > > +|7r+ > {tt- >) = — ^{V2\0,0 > +|2,0 >) (12.1) 
V 2 V 3 

J=odd 1^4 > = -^{{tt^ > |7r+ > -|7r+ > |7r~ >) = |1,0 > 
v2 

where we used the convention Ivr''" >= — 11,+1 > [57] and replaced the symmetrization normaliza- 
tion factor ^ by The states IS" > and \A > then acquire a meaning of maximally entangled 

Bell states of two-pion charge states with pion charge states | vr" > and vr^ > representing the two 
qubit states. 



The interaction with the environment is assumed to conserve the isospin Ik and G-parity. The 
G-parity of n pion states is G = (—1)" while isospin multiplet of particle-antiparticle pairs in LS 
states has G = (— 1)^+'S'+^ [53']. in vr'p —>■ ■K~TT~^n the qq and q^lf' resonant modes must be in 
triplet states of the quark and diquark spins [58] . The conservation of G-parity then implies that 
for resonant modes K + Ik = even. For states with J=even (odd) the transitions from i^=even 
(odd) will conserve Generalized Bose-Einstein symmetry while the transitions from i^=odd (even) 
will violate it. We thus can introduce Bose-Einstein symmetry conserving and violating two-pion 
isospin states 

|/c(vr-7r+) >= ^(Ivr- > \tt+ > +(-l)-^|7r+ > Ivr" >) (12.2) 



|/v(7r 7r+) >= -^{\tt > |7r+ > -(-l)'^|7r+ > |7r >) 
v2 

The summation in (11.6) over JX thus has two parts corresponding to two-pion isospin states 
Ic{Tr~TT~^) and Iy(7r~7r"^) and Generalized Bose-Einstein symmetry conservation and violation, 
respectively. 
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The two kinds of two-pion states isospin states (12.2) give rise to co-evolution amplitudes that 
conserve Generalized Bose-Einstein symmetry 

< JA(7r-7r+)r/,i"c(7r"7r+),r„| < 9u9d\T\gu9d > |0r >= (12.3) 
^ < JA(7r~7r+),/c(7r"7r+),r„| < gu9d\TfE\9u9d > \K n{qq) , Ic , tk ij. > 

axirr?) < KiJ,{qq)r], Ic,TKix\Ti\OT > 
and co-evolution amplitudes that violate Generalized Bose-Einstein symmetry 

< JA(7r"7r+)r?,/v(7r"7r+),rn| < 9u9d\T\9u9d > |0r >= (12.4) 

^ < JA(7r"7r+),/y(7r"7r+),r„| < 9u9d\TfE\9u9d > \K ii{qq) , ly , tr > 

aK{rr?) < Kn{qq)r],Iv,TK^\Ti\QT > 

where /c,/y=0 or 1 for |/c(7r~7r+) >, |Iy (7r~7r"'") >= \S > or 1^4 >, respectively. Since we focus 
on resonant modes we have omitted for the sake of brevity the non-resonant contributions with 
/if = 2 on r.h.s. of (12.3) for J=even and on r.h.s. of (12.4) for J=odd. 

The total co-evolution amplitude is a combination of the two sub-amplitudes (12.3) and (12.4) 
which in the most general form reads 

< JA(7r"7r+)ry,/;jA(7r"7r+),T„| < 9u9d\T\9u9d > |0t >= (12.5) 
ajx,ri < -''A(7r"7r+)?7,/c(7r"7r+),r„| < 9u9d\T\gu9d > |0t > -|- 

^J\,v < •^-^(7r"7r+)?7,/y(7r"7r+),r„| < 9u9d\T\gu9d > |0t > 

where the symbol £'jA(vr~7r+) refers to the fact that the observed two-pion states are in general 
entangled isospin states 

< JA(7r-7r+)r?, Ejx{tt-tt+)\ = (12.6) 

ajx,r, < JA(7r~7r+)77,/c(7r~7r+)| +ujjx,r, < JA(7r~7r+)?7, /y (7r~7r+)| 

The entanglement amplitudes 

ajx,v = < JA(7r-7r+)r?,/;jA(vr-7r+)|JA(7r-7r+)r?,/c(7r-7r+) > (12.7) 
i^JX,r, = < JA(7r"7r+)??,/;jA(7r"7r+)|JA(7r"7r+)?7,/y(7r"7r+) > 

are normalized 

\ajX,r,f + |wjA,r,P = 1 (12.8) 
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To see explicitely the entanglement of the two-pion charge states we rewrite (12.6) in the form 

< JA(7r-^+)r?, SjA(^~vr+)| =< JA(7r-7r+)r?| (ajA,^ < 7r-7r+)| + (-1)^6ja,„ < vr+7r-)|) (12.9) 
where 

ajX,V = -^{oiJX,r]+i^JX,r]) (12.10) 
bjX,rj = -^{otJX,ri - ^JX,r) 

For ajx^ri = ^jx,ri = the two-pion isospin state is separable. 

Nucleon helicity amplitudes with definite dipion spin J and helicity A are a combination of 
nucleon transversity amplitudes with the same spin J and helicity A but with opposite transversity 
T [i^, [4^. For any combination of solutions of transversity amplitudes the helicity amplitudes 
must share the same isospin state as the transversity amplitudes. The entanglement amplitudes 
thus cannot depend on transversity r and the quantum numbers gr- 

The plane wave co-evolution amplitudes describe the angular distribution of the two-pion states. 
Their angular expansion in terms of angular amplitudes reads [45] 

< 6'0,?7, S(7r"7r+),r„| < gu9g\T\gugd > |0r >= (12.11) 
"^¥^(9(1)) < JA?7,^(7r"7r+),r„| < gugg\T\gu9d > |0r >= 

J, A 

^Y^^{e(f>) < JXr],Ejx{7r^TT+),Tn\ < gu9g\T\gu9d > |0r > 

J, A 

where 6, 4> describe the direction of vr^ in the two-pion center-of-mass system. Assuming that the 
final two-pion isospin state \E{TT~ir^) > is entangled state of symmetric and antisymmetric states 

\E{TT-7r+)>=as\S > +aA\A> (12.12) 

where |a5p + la^P = 1, the self-consistency of the angular expansion (12.11) requires 

J=even: ajx,-q = as, wja,?? = cla (12.13) 
J=odd: ajA,,, = aa, wja,,, = as 

We call the entanglement of the state |£'(7r~7r~*") > dynamic because it is produced in a hadron 
interaction process. In general, the entanglement amplitudes 05 and may depend on energy s, 
dipion mass m and momentum transfer t. 

As the result of the interaction with the environment, there is a change of the entanglement 
content of the two-pion charge, or isospin, states with definite spin J from maximally entangled 
states IS* > or 1^ > required by the Generalized Bose-Einstein symmetry with separable 7r~7r^ 
isospin states in the final TT~Tr~^n state. The entanglement change comes entirely from the violation 
of Generalized Bose-Einstein symmetry and demonstrates the fact that the interaction of pion 
creation process with the environment is a non-dissipative dephasing interaction. Below ~ 1000 
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MeV where S- and P-wavea mplitudes dominate it is the environment induced p^{770) — /o(980) 
mixing that is responsible for the violations of Generalized Bose-Einstein symmetry and the 
entanglement of the ■7t~tt~^ isospin states in the final state. 

In a sequel paper we will show that the requirement of Generalized Bose-Eistein symmetry 
leads to three relations among partial wave intensities in 7r~7r"'", vr'^vr'^ and tt~^tt~^ production 
for even dipion spins. Available data violate these relations for 5- wave and D-wave intensi- 
ties providing a direct experimental evidence for a violation of Generalized Bose-Einstein symmetry. 

In standard Quantum Field Theory the initial and final states of particles are separable Fock 
states. CPT invariance requires that the amplitudes describing reaction 7r~p TT~TT~^n also 
describe the CPT conjugate reaction 7r"'"7r~n — > n'^p. This is possible since the initial state 
7r"^7r~n is separable and thus experimentally preparable. The underlining assumption that makes 
this statement of CPT invariance possible is that the two reactions are both isolated events in the 
Universe. As the result, the evolution of the pion creation process is unitary. 

The evidence for evolution from pure initial states to mixed final states in Ti~p 7r~7r+n 
presented in our previous work [45] necessitates to view the pion creation process as an open 



quantum system interacting with a quantum environment. According to Wald Theorem [34], this 
interaction must violate CPT symmetry. In this work we have reached the conclusion that the 
interaction with the environment leads to entanglement of tt~tt~^ charge, or isospin, states. The 
entanglement content carried by the produced two-pion states depends on the kinematics and 
dynamics of the pion creation process and the final states are no longer separable Fock states. 
The actual final states < Ocp, E{7r~TT^),Tn\ < Qugdl do not posses prepareable CPT conjugate 
states due to the entanglement of 7r~7r+ pairs and because the environment states IguSd > do not 
have well defined charge conjugate states. As a result the concept of CPT symmetry looses its 
meaning. 

The dynamic entanglement of final states is a distinct feature of CPT violation not limited 
to pion creation processes. Recently Bernabeu, Mavromatos and Sankar have shown that a CPT 
violating interaction of free maximally entangled neutral mesons M°lvr (e.g. K^K or B^B 
pairs) with an environment of fluctuations of space-time metric (space-time foam) will change the 
the entanglement of the M^M pairs and violate Generalized Bose-Einstein symmetry 4l|, |42] 

\E{M^jf) > = ^(|M° > |M° > -|M° > \jf >) (12.14) 
v2 

+ ^{\jf > \M^ > +\M^ > \lf >) 
v2 

The complex amplitude w is a CPT violating parameter associated with the term violating Gen- 
eralized Bose-Einstein symmetry. It is similar to our amplitude ujjx^r] in (12.6). However, to 
arises from CPT violating interaction of the propagating free M^'M^ pairs with environment over 
macroscopic distances while toj\^fj arises from CPT violating interaction of propagating resonant 
modes with environment inducing their transitions to two-pion states. Recent measurements of 
u by KLOE Collaboration [iJ] found small uj consistent with a zero. Nevertheless, the work of 
Bernabeu, Mavromatos and Sankar is the first attempt to explicitely relate entangled final states 
to quantum gravity. Their work opens the possibility that p°(770) - /o(980) mixing observed in 
the CERN measurements of tt~p — > TT~TT^n on polarized targets also arises from the low energy 
manifestations of quantum gravity. 



42 



XIII. CONCLUSIONS. 



The measured S- and P-wave density matrix elements form an autonomous subspace of 
reduced density matrix which is analytically solvable at any dipion mass m. The presence of 
p°(770) in the S'-wavc spectra jS'p arises from the p^{770) peak in the data component ai + 02 
which survives the subtraction of p^{770) peak in the P-wave amplitude 3|Lp in the relation 
|5p = ai + a2 — 3|i^p. The presence of /o(980) in the P-wave amplitudes l-L^P and |C/rfp arises 
from the structures at /o(980) mass in all other data components. The p°(770) — /°(980) mixing 
results from information on dynamics encoded in all measured density matrix elements which also 
encode the level splitting of the mass spectra. The analytical form of the level splitting of mass 
spectra reveals the existence of a new quantum number g characterizing the quantum states of 
the environment and allows to identify the four sets of solutions for the transversity amplitudes 
with the four co-evolution amplitudes required by the Kraus representation of final state density 
matrix in tt~p — > 7r~7r"*"n. 

A model of the CPT violating and non-dissipative interaction of the pion creation process 
with the environment proposes to explain the /9°(770) - /o(980) mixing and predicts a dynamic 
entanglement of 7r~7r+ pairs in Tr~p 7r~7r + n. Theoretical work by Wald, Hawking, Ellis, 
Mavromatos and others suggests that the origin of the environment and its CPT violating non- 
local interactions is in quantum gravity rendered observable at low energies by the pion creation 
process. 
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